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ABSTRACT 

A  significant  part  of  radar  technology  has  been  concerned 
with  developing  techniques  which  permit  satisfactory  opera- 
tion of  radar  systems  when  moving  targets  are  in  the  presence 
of  relatively  fixed  targets  (clutter) , 

This  thesis  presents  the  significant  aspects  of  old, 
improved  and  modern  anticlutter  techniques,  each  of  which  is 
described  with  its  principal  advantages  and  disadvantages. 

A  tabulation  is  made  which  compares  each  anticlutter 
technique  on  the  basis  of  its  clutter  suppression  capability. 

Continuing,  an  application  of  digital  filters  is  made  by 
designing  a  fourth-order  recursive  filter  with  Butterworth 
characteristics,  using  a  computer  program.   This  filter  can 
be  used  in  a  range-gated  MTI  (Moving  Target  Indication) 
receiver  to  improve  its  performance. 

Special  consideration  was  given  to  the  AN/UPS-1  and 
AN/SPS-40B,C  radars  which  are  used  in  the  Greek  Navy  in  a 
clutter  environment  characteristic  of  the  Mediterranean 
Sea.   Appropriate  modifications  to  their  existing  analog 
systems  are  suggested. 
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I.   INTRODUCTION 

The  conglomeration  of  unwanted  echoes  coming  from 
fixed  targets  is  defined  as  "clutter".   These  echoes  clutter 
the  radar  display  and  make  the  recognition  of  wanted  signals 
difficult.   The  result  degrades  the  quality  of  radar 
information. 

There  are  several  types  of  clutter  such  as  those  caused 
by  weather,  ground,  sea,  birds,  chaffs,  rain,  clouds,  preci- 
pitation and  other  meteorological  phenomena.   In  addition, 
second-time-around  clutter  is  often  troublesome. 

For  these  types  of  clutter,  there  are  a  number  of  reme- 
dies for  improving  the  signal-to-clutter  ratio.   Consequently 
radar  systems  are  configured  using  one  or  more  of  these 
methods. 

The  basic  signal  parameter  which  is  used  for  distin- 
guishing fixed  and  moving  targets  is  the  doppler  shift  in 
frequency  of  the  moving  target.   This  frequency  shift  is  given 
by  the  relationship 

2V 

f^  =  -^  (1-1) 

where  V   is  the  radial  relative  velocity  component  of  the 
target  and  X    is  the  wavelength  of  the  transmitted  signal. 
The  presence  of  various  frequencies  is  the  basis  for 
the  use  of  a  filter  that  passes  echoes  with  a  doppler  shift 
but  rejects  clutter  which  has  zero  or  small  doppler  shift. 


Over  the  last  twenty  years,  a  concerted  effort  has  been 
made  to  understand  the  root  causes  of  clutter,  its  charac- 
teristics, and  to  find  ways  for  MTI  performance  to  be 
improved.   This  effort  has  resulted  in  new  signal  processing 
techniques,  which  together  with  certain  antenna  arrange- 
ments, minimize  clutter  problems.         ^. 

Analog  MTI  is  one  technique  which  has  been  used  in  radar 
systems  for  a  number  of  years  by  the  Greek  Navy  with 
limited  success.   Performance  has  suffered  from  the  inherent 
difficulty  of  controlling  the  delay  time  and  system  gains. 
The  on-line  availability  has  been  generally  poor.   The  radars 
operate  well  when  tuned-up,  but  experience  deteriorating 
performance  with  time.   The  effective  contribution  to  the 
overall  radar  system  performance  has  generally  been  limited 
due  to  the  inflexibility  of  the  hardware. 

However,  many  of  the  disadvantages  of  the  analog  MTI 
have  been  eliminated  by  digital  MTI.   This  is  shown  in  this 
thesis  using  computer  programs  to  design  an  optimum  filter. 
No  attempt  was  made  to  implement  the  filter  with  hardware. 
However,  the  computer  results  are  in  a  form  consistent  with 
the  format  needed  for  easy  digital  filter  design. 

It  is  expected  that  when  the  entire  system  is  built  and 
operated,  it  will  have  the  inherent  flexibility  and  stability 
required  to  obtain  greatly  improved  performance  against 
clutter . 


II.   OLD  TECHNIQUES 

In  the  design  of  the  original  radar  systems,  the  main 
effort  was  directed  toward  the  ability  of  the  radar  receiver 
to  detect  the  presence  of  echoes  and  extract  information  from 
them. 

Fundamentally,  detection  was  limited  by  the  presence  of 
noise  in  the  receiver.   In  addition  detection  of  moving 
targets  was  inhibited  by  the  presence  of  clutter  frequently 
much  stronger  than  the  moving  targets  themselves. 

Under  such  circumstances  the  radar  became  useless.   The 
need  to  overcome  this  led  to  the  design  of  additional  analog 
subsystems  (techniques) ,  which  for  that  time,  offered 
considerable  detection  improvement.   A  number  of  these  could 
be  mentioned,  but  those  few  which  are  still  employed  in  older 
but  existing  radar  systems  will  be  described. 

Part  II  covers  the  following  techniques. 

(i)  Frequency  Agility 

(ii)  Circular  Polcarization 

(iii)  Sensitivity  Time  Control  -  Logarithmic  Receivers 

(iv)  Instantaneous  Automatic  Gain  Control 

A  brief  description  of  each  is  given,  indicating  the 

advantages  and  disadvantages. 

A.   FREQUENCY  AGILITY 

Frequency  agility  is  a  technique  which  enables  one  to 
achieve  a  measure  of  clutter  reduction  by  changing  the 


frequency.   Coaxial  magnetrons  permit  a  change  in  the  basic 
frequency  by  mechanical  means.   Signal  enhancement  in  clutter 
by  a  frequency  change  is  shown  in  Figure  1  which  is  an  A-scope 
presentation  of  a  rooftop  radar  [Refs.  1/2]. 

Part  (a)  of  Figure  1  [Ref.  1]  is  with  the  agility  switch 
ON  and  part  (b)  OFF.   The  radar  was  operated  to  look  at  small 
targets  over  a  fairly  cluttered  ground  terrain.   In  this 
case  the  bandwidth  of  agility  was  about  75  MHz.   Clearly  a 
distinct  improvement  is  seen  and  the  two  hard-targets  are 
more  easily  discerned. 

If  there  is  a  doubt  about  the  presence  of  targets,  then 
a  quick  flip  of  the  Agility  switch  from  ON  to  OFF  can 
usually  help  resolve  the  doubt. 

To  obtain  frequency  agility  mechanical  adjustments  are 
made  within  the  microwave  tube.   These  adjustments  may  be 
as  follows: 

(i)    By  moving  the  end  wall  of  the  coaxial  cavity  up 

and  down  slowly.   This  movement  results  in  broad- 
band tuning, 
(ii)   By  varying  the  ring's  diameter  of  the  outer 

periphery  of  the  cavity.   This  form  is  independent 
of  (i)  and  causes  changes  in  the  tuning  sensitivity 
of  the  cavity. 
The  second  technique  gives  a  tuning  speed  ranging  from  80 
to  200  Hz/second  and  a  tuning  bandwith  up  to  100  MHz. 
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Figure  1.   A-scope  presentations,  (a)  with  and 
(b)  without  frequency  agility. 
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The  advantage  of  this  type  of  tuning  is  the  improved 
performance  as  demonstrated  by  its  use  in  operational  systems 
today. 

However,  there  is  a  disadvantage,  when  precipitation 
extends  over  a  large  volume  and  completely  fills  the  antenna 
beam.   Then,  the  total  energy  received  from  all  the  preci- 
pitation particles  within  the  radar  resolution  cell  can  be 
quite  large  compared  with  that  received  from  a  conventional 
target  and  can  mask  or  clutter  the  radar  display. 

In  this  case  the  amount  of  clutter  with  which  the 
target  echo  must  compete  may  be  reduced  with  narrow  beam- 
widths  and  narrower  pulse  widths  if  otherwise  acceptable. 

The  best  method,  which  can  be  applied  to  any  radar 
system,  of  eliminating  weather  effects,  such  as  heavy 
precipitation,  is  to  operate  at  lower  radar  frequencies. 
It  may  be  noted  that  in  weather  radar  applications,  the 
designer  has  the  option  of  operating  at  the  lower  frequen- 
cies which  give  a  better  S/C  ratio. 

B.   CIRCULAR  POLARIZATION 

When  a  radar  operates  at  the  higher  frequencies  and 
the  target  is  of  a  complex  nature  like  an  aircraft,  it  is 
possible  to  improve  the  detection  capability  in  the  presence 
of  precipitation  by  radiating  circularly  polarized  energy 
[Ref.  3]. 

Successful  detection  is  based  on  the  dissimilarity 
between  clutter  and  aircraft  echoes.   With  precipitation 
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and  raindrops,  the  backscattered  energy  does  not  depend 
on  the  polarization  of  the  incident  waves.   But  in  the  case 
of  the  aircraft  target  the  energy  is  dependent  on  the 
incident  polarization.   For  example,  the  radar  may  use  an 
antenna  which  is  responsive  to  rotation  of  the  polarized 
wave  in  one  direction.   Since  symmetrical  reflectors  reverse 
the  rotation  of  the  wave,  that  echo  energy  will  not  be 
accepted  by  the  receiver. 

Circular  polarization  is  generated  by  two  methods  as 
outlined  in  Ref.  3  and  as  shown  in  Figure  2(a,b).   The  first 
part  shows  a  quarter-wave  plate  which  converts  linear  polari- 
zation to  circular  by  passing  the  wave  through  it.   The 
second  part  shows  a  turnstile  junction  which  is  a  combination 
of  rectangular  and  circular  waveguides  properly  matched  to 
generate  circular  polarized  waves. 

Generally  with  circular  polarization  the  ratio  of  target 
to  precipitation  signal  is  improved  by  8  to  25  dB  over  that 
obtained  with  linear  polarization.   Under  normal  conditions 
an  improvement  of  15  to  20  dB  is  expected.   It  must  be  noted 
that  in  all  previous  cases  the  ground  reflections  are 
considered  negligible. 

C.   SENSITIVITY  TIME  CONTROL  -  LOGARITHMIC 
RECEIVER  CHARACTERISTICS 

1.   Sensitivity  Time  Control  (STC) 

One  of  the  earliest  and  most  widespread  but  less 

sophisticated  techniques  for  attenuating  the  clutter  signal 
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Figure  2.   Quarter-wave  plate  and  turnstile  junction 

for  converting  linear  to  circular  polarization 
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is  the  gain  at  close  ranges.   This  is  known  as  STC  [Refs.  2, 
1,    13] .   The  reduction  of  clutter  from  close  targets  is 
achieved  by  reducing  the  likelihood  of  receiver  saturation, 
which  leads  to  a  programmable  receiver  gain  as  a  function  of 
time  so  that  the  gain  will  be  low  for  targets  at  close 
range  and  high  for  targets  at  long  range. 

The  problem  of  dynamic  range  limitation  of  radar 
receivers  was  overcome  by  switching  on  the  STC  circuit. 
This  simply  reduced  the  gain  of  the  amplifiers,  usually  at 
the  intermediate  frequency  stage.   The  gain  varies  as  R  , 
where  n  can  take  a  value  between  2  for  weather  effects  and 
4  for  land  and  sea  clutter. 

The  resulting  improvement  appears  on  the  PPI,  where 
the  STC  circuit  reduces  the  excessive  brightness  at  the 
center  of  the  display  due  to  the  bunching  of  the  radial 
sweep  lines.   This  anti-clutter  circuit  is  used  only  where 
needed  and  should  be  switched  off  at  other  times.   Conse- 
quently, the  use  of  this  technique  depends  upon  an  operator's 
judgment  and  is  more  helpful  the  more  experienced  the 
operator. 

In  addition  to  the  above,  this  type  of  anti-clutter 
technique  has  the  limitations  caused  by  its  analog  circuitry 
and  IF  stage  processing,  which  limit  the  radar's  clutter 
reduction  and  target  detection  capabilities. 

Similar  to  this  is  the  fast-time  constant  (FTC) , 
another  analog  circuit.   This  is  normally  a  high-pass 
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filter  (RC) ,  as  outlined  in  Ref.  3.   Its  use  also  depends 
on  the  judgement  of  the  radar  operator.   It  is  coiranonly 
used  in  conjunction  with  STC  and  is  effective  on  all 
targets  including  clutter  at  all  ranges  of  the  radar. 
2.   Logarithmic  Receiver  Characteristics 

The  receiver  gain,  when  clutter  echoes  are  present, 
can  be  reduced  to  prevent  relatively  large  clutter  signals 
from  saturating  the  receiver  by  using  a  logarithmic 
amplifier  followed  by  a  differentiating  circuit  or  a 
pulse-length  discriminator  [Refs.  2,    3]. 

This  technique  is  used  to  reduce  the  effective 
dynamic  range  of  clutter  or  jamming  signals  at  the  threshold 
detector  or  display.   Since  clutter  varies  in  intensity  with 
range,  the  use  of  a  log-receiver  reduces  this  and  provides 
a  clutter-signal  amplitude  independent  of  range  at  the 
output. 

The  transfer  characteristic  of  a  log-amplifier  has 
the  form 

y   =   a  ln(l  +  bx)  (2.1) 

where  a,  b  are  constants  of  the  amplifier,  x  is  the  instan- 
taneous amplitude  of  the  input  and  y  is  the  output  [Ref.  2J . 

The  purpose  of  the  differentiator  circuit,  which 
follows  the  log-amplifier,  is  to  reduce  the  mean  value  of 
the  noise  return  or  clutter.   This  is  successful  only  when 
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the  clutter  has  a  uniform  distribution  in  range  for  a  period 
greater  than  that  of  the  FTC  circuit. 

At  this  point,  techniques  1  and  2  are  compared. 
Sensitivity  Time  Control  (STC)  suppresses  clutter  echoes, 
which  enter  via  the  sidelobes  much  better.   STC  automatically 
turns  down  the  gain  at  close  ranges,  reducing  in  that  way 
sidelobe  signals.   The  log-amplifier  does  not  suppress  side- 
lobe  clutter  which  masks  the  radar  display.   Both  techniques 
have  the  disadvantage  of  requiring  perfect  adjustment  and 
skilled  operators. 

There  are  a  few  applications  where  these  two  tech- 
niques may  be  employed  together.   Marine  radars  use  log- 
receivers  exclusively  because  they  permit  strong  and  weak 
echoes  (returns  of  sea,  buoys  and  boats)  to  appear  on  the 
radar  scope. 

The  block  diagram  of  Figure  3  corresponds  to  a  log- 
amp  followed  by  a  pulse-length  discriminator  (PLD) .   This 
form  of  receiver  is  an  alternative  solution  having  the  same 
effect  on  clutter  and  noise  as  the  log-FTC  circuitry, 
rejecting,  as  well,  forms  of  land  clutter  and  pulsed 
interference . 

The  bandpass  filter  plays  the  role  of  the  matched 
filter  of  the  conventional  linear  receiver,  and  passes 
pulses  within  fixed  frequency  limits  which  are  equal  to  the 
pulse  period  of  the  delay  line. 

The  disadvantage  of  this  technique  is  the  same  as 
the  log-FTC,  in  that  loss  of  detectability  results.   These 
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losses  reach  the  level  of  4-8  dB  for  single-hit  detection 
and  2-4  dB  for  u-Pulse  postdetection  integrations. 

D.   INSTANTANEOUS  AUTOMATIC  GAIN  CONTROL  (lAGC) 

The  lAGC  technique  is  a  fast-acting  form  of  AGC  to 
prevent  the  saturation  of  IF  amplifiers.   It  is  used  in 
most  radar  receivers  and  acts  in  a  similar  manner  to  the 
automatic  volume  control  (AVC)  in  a  radio  receiver. 

A  block  diagram  is  shown  in  Figure  4,    with  greater 
detail  given  in  Ref .  3.   The  control  of  the  gain  of  the  IF 
amplifier  is  obtained  by  using  negative  feedback.   Its 
filter  has  a  narrow  bandwidth  to  minimize  the  effect  of  the 
extended  clutter  echoes. 

There  is  a  limitation  on  the  lAGC  loop,  which  can  only 
control  one  or  two  stages.   For  this  reason  instabilities 
sometimes  appear  in  a  system  with  more  stages. 

A  second  technique  similar  to  lAGC  is  that  known  as 
detector  balanced  bias  (DDB) ,  which  is  shown  in  the  diagram 
of  Figure  5.   DDB  differs  from  lAGC  in  three  basic  respects; 

(i)      Does  not  employ  a  feedback  loop. 

(ii)     Bias  is  applied  to  a  diode  detector  instead 
of  the  IF  amplifier. 

(iii)    Employs  a  delay  line  in  the  absence  of  the 
feedback  loop. 

While  some  of  these  techniques  used  singularly  or  in 
combination  are  beneficial  in  reducing  the  effects  of  some 
types  of  clutter,  they  are  not  sufficient  in  modern 
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Figure  5.   Detector  balanced  bias  (DDB) 
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sophisticated  radar  systems.   More  importantly  they  are  not 
effective  in  reducing  the  types  of  clutter  that  occur  in 
the  Mediterranean  environment  which  is  the  concern  of  this 
thesis . 
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III.   IMPROVED  TECHNIQUES 

In  World  War  II  a  particular  class  of  radar  processors 
made  its  appearance.   This  type  of  processor  extracted  from 
the  echo  the  doppler  shift  in  frequency  caused  by  a  moving 
target.   Originally,  such  radar  systems  employed  CW  trans- 
mission.  The  newer  systems  were  able  to  discriminate 
between  the  echoes  of  moving  and  fixed  targets. 

At  that  time  the  problem  of  clutter  suppression  was 
approached  from  many  directions;  "clutter"  meaning  the 
return  echoes  from  fixed  targets  in  the  vicinity  of  the 
target  of  interest.   Sea,  ground,  weather  effects,  angels, 
birds,  rain,  chaff,  trees  and  vegetation  movement  are  common 
sources  of  clutter. 

The  analog  techniques  originally  developed  usually 
employed  delay  line  cancelers  or  analog  filter  banks  (matched) 
These  techniques  are  still  used  in  large  numbers  in  contem- 
porary applications.   Delay  line  cancelers  are  used  in  the 
UPS-1  and  AN/SPS-40B,C  radars  on  ships  of  the  Greek  Navy. 

Thus,  this  thesis  continues  to  compare  the  performance 
of  other  techniques,  both  existing  and  modern,  in  an  effort 
to  find  a  satisfactory  clutter  rejection  filter  that  is  not 
unduly  complex  and  can  be  rather  easily  implemented. 

Part  III  describes  the  following  representative  tech- 
niques with  their  respective  advantages  and  disadvantages. 

(i)      Analog  MTI  and  Pulse-Doppler  Systems. 
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(ii)     Dual  Beam  Antenna, 
(iii)    Log  -  CFAR. 
(iv)     Video  Detector. 

A.   ANALOG  MIT  AND  PULSE-DOPPLER  SYSTEMS 

When  doppler  information  is  extracted  from  a  pulsed 
radar  return,  it  is  known  as  an  MTI  radar.   The  MTI  stands 
for  "moving-target  indication".   It  may  also  be  called 
pulse-doppler  radar.   In  practice,  a  distinction  is  some- 
times made  between  the  MTI  radar  and  the  pulse-doppler  radar, 
although  they  are  both  based  on  the  same  principle. 

MTI  usually  refers  to  a  radar  in  which  the  doppler 
frequency  measurement  is  ambiguous  but  the  range  measurement 
is  unambiguous.   Another  characteristic  feature  of  the  MTI 
radar  is  the  delay-line  canceler  used  to  detect  the  doppler 
frequency  shift. 

In  the  pulse-doppler  radar  the  doppler  measurement  is 
usually  unambiguous  and  the  range  may  or  may  not  be  ambiguous 
Ambiguous  range  means  that  multiple-time-around  echoes  are 
possible,  while  ambiguous  doppler  implies  that  "blind  speeds" 
fall  within  the  range  of  expected  target  speeds. 

MTI  radar  has  the  ability  of  extracting  the  moving-target 
echo  from  the  clutter  echo  even  when  the  clutter  is  20  to 
30  dB  greater  than  the  moving-target  echo.   Some  pulse 
doppler  radars  can  detect  moving  targets  even  when  the 
clutter  echo  is  70  to  90  dB  greater  than  the  target  echo. 
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The  distinction  between  the  two  radars  arose  histori- 
cally, and  is  commonly  still  made.  In  many  instances  the 
difference  between  them  is  only  a  matter  of  nomenclature. 

Figure  6  shows  the  simplest  MTI  processor,  the  single 
delay-line  canceler,  which  extracts  the  doppler  information 
by  subtracting  two  successive  echoes.   The  delay-line 
canceler  plays  the  role  of  a  filter,  eliminating  the  d-c 
component  of  fixed  targets  and  permitting  the  a-c  components 
of  moving  targets  to  pass. 

There  are  several  types  of  MTI  processors.   A  distinc- 
tion between  them  is  the  type  of  information  processed  in 
the  returned  signals,  that  is,  whether  the  phase,  the 
amplitude,  or  both  phase  and  amplitude  are  processed.   These 
types  are  described  with  their  block  diagrams  in  Refs.  2  and 
3.   It  is  beyond  the  stated  intention  of  this  thesis  to 
describe  all  types  of  MTI,  but  some  are  listed  here  for 
reasons  of  interest. 

(i)     Phase-Processing  MTI 

(ii)    Amplitude-Processing  MTI 

(iii)   Vector-Processing 

(iv)    Clutter-Locking  MTI 

(v)     Multiple-Canceler  Systems 

These  processors  do  not  satisfactorily  handle  the  well 
known  problems  of  clutter  precipitation.   This  problem 
becomes  more  important  when  the  effects  of  wind  shear  and 
average  clutter  velocity  are  taken  into  account.   For  example. 
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the  clutter  performance  of  the  aforementioned  three  radars 
is  poor  under  these  conditions.   This  has  been  verified  by 
the  experience  of  the  Greek  Navy  in  the  Mediterranean 
environment . 

To  get  a  clearer  idea  of  this  limitation  on  MTI  per- 
formance, the  clutter  improvement  factor  I  was  plotted  vs. 
range.   Figure  7  shows  the  plot.   The  equation  of  I  is 
as  follows 

2   2 

8   (a:^u^j3  +  .18k^R%^) 

and  the  Continuous  System  Modeling  Program  1  was  used  with 
some  typical  values  given: 

a,   ,   =   variance  of  wind  turbulence,  varying  from 

.7  to  1.1  meter/sec 

X     =   .30  m,  transmitted  wavelength 

f      =   600  Hz,  PRF 
r 

k     =5.7  meter/ (sec) (km),  is  the  component  of  wind 
velocity  gradient  in  the  direction  of  the 
radar  beam 
R     =   range  in  km 

(ft     =   1.5°,  vertical  half -power  beamwidth-two-way 
path  in  radians. 
Figure  7  shows  that  the  attenuation  of  clutter  is  quite 
good  for  ranges  up  to  18  km  where,  for  <?    j^  =  .7  and 
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Figure  7.   Clutter  improvement  factor  vs.  Range 
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1.1  meter/sec  a  difference  of  1.5  dB  is  measured.   As  range 
increases,  the  clutter  attenuation  decreases  rapidly  as 
indicated. 

B.   DUAL  BEAM  ANTENNA 

The  MTI  visibility  is  highly  degraded  at  short  ranges 
due  to  nearby  graound  clutter,  sea  clutter  and  angel  echoes. 
To  eliminate  these  problems  a  dual  beam  antenna  system  has 
been  designed. 

This  system  is  shown  in  Figure  8.   Using  the  low  and 
high  beam  horn  at  the  appropriate  time,  as  outlined  in 
Ref.  13,  the  clutter  is  reduced  at  the  receiver  input.   The 
change  in  the  use  of  horns  occurs  by  switching.   The  two 
horn  noses  have  a  difference  of  about  4  degrees  in  elevation, 
The  angle  is  variable. 

The  degree  of  clutter  suppression  and  target  visibility 
can  be  determined  by  the  beam  angle  and  switching  distance. 
The  coverage  is  varied  by  varying  the  angle  difference  in 
the  nose  directions.   Consequently,  the  switching  distance 
is  changed  and  this  determines  the  degree  of  clutter 
suppression. 

Figure  9  [Ref.  13]  compares  the  S/C  ratio  of  dual  beam 
antennae  with  single  beam  antennae.   After  determining 
target/clutter  ratios  by  a  computer  simulation  program/ 
conclusions  from  these  two  techniques  are  obtained.   This 
comparison  shows  that  a  better  improvement  factor  is 
obtained  using  a  dual  beam  antenna  system. 
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Figure    9.      Comparison   in  S/C   ratios   between  dual   beam 
and   single   beam  antenna   systems. 


30 


Types  of  clutter,  which  are  taken  into  account  here 
are  ground,  sea  and  angel  echoes.   In  all  cases  an  improve- 
ment of  15  to  20  dB  is  obtained  using  dual  beam  antennas. 
Detection,  however,  becomes  difficult  when  the  range  is  less 
than  5  nautical  miles. 

These  advantages  of  the  dual  beam  antenna  are  the  reasons 
which  make  it  useful  in  systems  such  as  ARSR  (Air  Route 
Surveillance  Radar) . 

C.   LOG-CONSTANT  FALSE  ALARM  RATE 

Additionally,  a  special  technique  for  ARSR  systems  is 
that  of  Log-CFAR,  which  has  been  widely  used  to  eliminate 
weather  clutter.   The  reason  this  technique  has  been 
employed  is  that  neither  MTI  nor  dual  beam  antenna  systems 
have  the  ability  to  suppress  weather  clutter  which  has  a 
radial  velocity  distribution  at  high  angles  of  elevation. 

Previously  mentioned  methods  obtain  some  improvement 
by  improving  the  signal-to-clutter  ratio.   But  in  the  case 
of  clutter,  whose  signal  strength  is  well  above  the  threshold 
or  noise  level,  the  output  of  the  receiver  is  saturated. 

An  additional  suppression  of  clutter  is  needed  and  is 
obtained  by  using  the  Log-CFAR  method.   It  keeps  a  constant 
output  at  the  receiver  suppressing  the  clutter  to  noise 
level . 

In  terms  of  the  block  diagram  this  technique  is  shown 
in  Figure  10  and  its  operation  is  outlined  in  Ref.  13.   A 
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step  by  step  explanation  is  given  in  terms  of  input-output 
waveforms  according  to  the  diagram  of  Figure  10. 
The  output  w  of  log-amplifier  is  given  by 

w   =   a  ln(l  +  bv)  (3.8) 

where  a,b  are  constants  and  v  is  the  input  [II-C-2] .   The 
subtracter  has  an  output  y,  which  is  given  by 


y   =  w  -  X   =   a  In  (^)  +  ^         (3.9) 


where  x  is  the  average  detector  value  given  by 


X 


=   ^   I   ^i  (3.10) 


and  where  y  is  the  Euler's  constant,  a  the  variance  of  the 
input  signal  v  and  N,  the  sample  number.   The  output  signal 
of  the  anti-log  converter  is  given  by 


=   c  exp  (d  .  y)   =   ^  exp  (^)      (3.11) 


where  c,d  are  constants  and  a  •  d  =  1. 

Sample  number,  N,  plays  an  important  role,  since  it 
determines  the  processing  loss  of  the  CFAR  circuit  and 
consequently  the  clutter  suppression.   Reference  13  shows 
that  the  cutoff  frequency  of  the  rejection  band  is  given  by 
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fc   =   JTT  (3-12) 


where  t  is  the  pulse  width. 

The  Log-CFAR  process  has  two  disadvantages  which  are 
listed  below. 

(i)     Sample  number  N  must  have  a  value  of  about  8 

to  keep  the  processing  loss  between  1  and  2  dB. 
(ii)    Mean  value  of  clutter  component  must  have  a 
frequency  less  than  l/Nx  to  be  eliminated. 

D.   VIDEO  DETECTOR 

This  application  is  specially  useful  in  non-coherent 
airborne  radars  operating  either  over  land  or  in  the  proximity 
of  a  land  mass,  where  very  intense  ground  clutter  patches  or 
individual  scattering  objects  exist. 

Such  ground  clutter  may  have  radar  cross-sections  on  the 
order  of  one  to  ten  million  square  meters.   Thus,  even  with 
the  best  designed  antennas  on  airborne  platforms,  this 
clutter  may  be  strong  enough  to  be  detectable  through  the 
antenna  sidelobes. 

These  signals  in  general  will  have  doppler  velocities 
within  the  acceptance  band  of  the  MTI  filter  or  filter  bank 
and  will  result  in  a  high  rate  of  false  alarms.   The  side- 
lobe  clutter  problem  is  a  very  serious  one  but  the  use  of 
a  square-law  detector  provides  significantly  improved  clutter 
rejection. 
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For  main  lobe  clutter  which  is  assumed  to  be  Gaussian 
in  both  amplitude  and  frequency,  the  MTI  improvement  factor 
is  calculated  for  different  spectral  widths  using  linear 
and  square  law  detector  characteristics. 

The  block  diagram  of  the  system  is  shown  in  Figure  11 
where  a  magnetron  is  used  as  the  transmitter.   The  echo 
signal  is  converted  to  IF  and  then  amplified.   Detection  is 
realized  by  the  so  called  Video  Detector.   After  MTI 
filtering,  digital  signal  processing  follows  before  the 
video  presentation  on  a  display. 

The  signal  x(t)  has  a  general  spectrum  in  the  IF  shown 
in  Figure  12.   Assuming  that  the  clutter  component  is  much 
larger  than  the  other  components ,    the  video  spectrum  is 
shown  in  Figure  13.   After  the  video  detector  this  spectrum 
is  approximately  the  auto-convolution  of  the  input  spectrum. 
Convolutions  not  containing  C,  which  stands  for  clutter,  are 
neglected  from  the  spectrum.   Component  C*C  is  undesirable 
clutter  and  is  suppressed  by  a  digital  high-pass  filter. 

The  system  described  above,  using  the  square  law 
detector  provides  a  wider  region  of  good  clutter  suppression- 
It  is  clearly  superior  in  comparison  with  other  improved 
clutter  suppression  techniques,  but  the  dynamic  range  of  the 
output  signal  is,  of  course,  larger  for  the  square  law 
detector  which  is  a  disadvantage  in  several  applications. 
It  can  be  noted  that  at  ranges  where  there  is  no  clutter 
the  target  is  not  detectable  after  the  MTI  filter.   For 
target  detection,  the  filter  must  be  bypassed  at  such  ranges. 
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II 

Clutter  beyond  the  unambiguous  range,  after  convolution  with 
first  time  around  echoes,  results  in  white  noise.   This  is 
an  additional  system  problem. 
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IV.   MODERN  TECHNIQUES 

This  section  is  prepared  to  demonstrate  several 
"modern  techniques"  which  are  employed  in  radar  system 
technology.   The  word  "modern"  covers  the  developments  of 
the  past  5  years  and  includes  future  trends  in  signal  wave- 
form design  and  processing  technology.   It  describes  three 
techniques,  which  are  listed  below. 

(i)      Pulse  Doppler  Processing 

(ii)     Moving  Target  Detector 

(iii)    Application  of  FFT  Algorithm. 

Finally,  a  tabulation  is  made  to  compare  the  amount  of 
clutter  suppression  achieved  by  each  technique  in  the  three 
categories  on  different  types  of  clutter. 

A.   PULSE-DOPPLER  PROCESSING 

Here  a  technique  is  presented  which  is  especially 
applicable  to  high-power,  pulse-doppler  radars,  achieving 
good  clutter  discrimination.   The  basis  of  this  method  is 
an  adaptable  filtering  and  integration  system  which  processes 
the  signal  in  a  way  that  the  processin   loss  is  minimized. 
Besides,  good  range  resolution  and  velocity  resolution  can 
be  obtained.   The  integration  is  made  adapative  to  the 
changes  in  the  variable  PRF  by  using  commutated- capacitor 
filters. 
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The  clutter  filter  system  has  a  block  diagram  shown 
in  Figure  14.   The  signal  waveform  looks  like  a  Gaussian- 
modulated  envelope  due  to  its  antenna  characteristics.   A 
Gaussian  filter  is  placed  at  the  beginning,  which  broadens 
the  duration,  t,  of  the  pulses  and  consequently  reduces 
their  peak  amplitude. 

This  filter  provides  a  narrow  spectriim  of  the  signal 

for  adjacent  channel  rejection  during  the  mixing  processes. 

Then  a  Gaussian  low  pass  filter  rejects  the  fundamental 

frequency  f   and  all  harmonics  generated  by  the  previous 

process.   After  the  sample  and  hold  circuit  a  high  pass 

filter  is  connected.   The  purpose  of  this  filter  is  to 

attenuate  clutter  components  to  below  noise  level. 

The  two  signals  are  then  mixed  with  f  .   Since  these 
^  o 

are  out  of  phase  by  tt/2  radians,  they  are  combined  to  provide 
an  output  signal  from  which  clutter  has  been  removed. 

Some  important  advantages  of  this  technique  are  listed 
below  [Ref .  16] . 

(i)       Clutter  suppression  at  levels  of  about  80  dB 

down  in  amplitude. 
(ii)      Use  of  capacitors  in  association  with 

adaptable  sample  and  hold  circuits,  stationary 
clutter  is  reduced, 
(iii)     Rapid  changes  of  PRFs  during  the  time  on 

target  permits  the  system  to  respond  without 
spreading  the  clutter  or  destroying  the 
integrity  of  the  signal  information. 
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(iv)      The  low-pass  filter  reduces  system  noise 
to  a  much  lower  level. 

The  disadvantages  of  the  system  are  the  range  and 
velocity  ambiguities  which  occur  when  it  is  required  to 
detect  very  fast  targets.   Also,  the  greater  complexity 
due  to  three  distinct  stages  of  spectrum  analysis. 

This  study  leads  to  the  conclusion  that  the  introduction 
of  clutter  filters,  which  are  switched  in  synchronism  with 
the  pulse  repetition  frequency  sequence,  effectively  reduce 
the  generation  of  clutter  side-bands  by  an  additional 
40-50  dB  over  a  fixed  filter  system.   It  is  thus  possible 
to  design  a  long  range  Pulse  Doppler  Radar  System  with  more 
than  80  dB  clutter  suppression,  if  the  PRFs  can  be  switched 
at  a  rapid  rate . 

The  storage  elements  are  capacitors  which  are  switched 
by  a  hardware  program  to  integrate  or  discharge  signals. 
By  this  means  real-time  target  assessment  is  realized.   The 
performance  of  this  radar  system  is  consequently  significantly 
improved. 

B.   MOVING  TARGET  DETECTOR 

Recently,  the  Lincoln  Laboratory,  associated  with 
Massachusetts  Institute  of  Technology,  developed  a  signal 
processor  named  the  Fast  Digital  Processor  which  performs  a 
Digital  Fourier  Transform  on  the  received  signal  and  presents 
the  Fourier  coefficients  [Refs.  1,  17,  18]. 
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This  processor  is  very  fast,  since  it  is  a  special- 
purpose  computer  designed  to  perform  signal  processing  tasks 
such  as  FFTs.   Its  algorithm  was  tested  by  embodying  it  in 
a  hardwired  version  of  the  Moving  Target  Detector  (MTD) . 
The  MTD  has  been  designed,  built  and  tested  at  Lincoln 
Laboratory  where  it  is  now  undergoing   further  tests. 

The  MTD  technique  uses  the  cascade  of  an  MTI  and  pulse 
doppler  processor  to  achieve  clutter  rejection.   This 
combination  is  an  approximation  to  an  optimum  clutter 
rejection  filter. 

In  terms  of  the  block  diagram,  the  MTD  is  shown  in 
Figure  15,  where  it  is  seen  that  the  signal  to  be  processed 
arrives  at  the  IF  stage  from  the  output  of  the  IF  pre- 
amplifier and  is  fed  through  a  special  wide-dynamic-range 
amplifier  to  the  quadrature  video  detectors. 

The  two  quadrature  video  detector  outputs  are  converted 
to  10-bit  digital  numbers  by  the  A/D  converters.   The  system 
contains  an  8000-word  input  memory  and  about  900  integrated 
circuits.   A  disc  memory  is  used  as  a  fine-grained  ground 
clutter  map.   This  memory  allows  the  removal  of  such 
returns  by  the  use  of  adaptive  thresholds  in  range,  doppler 
and  angle. 

It  is  required  that  a  number  of  scan  returns  be  observed 
before  a  target  is  declared  present  and  track  is  initiated. 
Thus,  track  initiation  is  used  to  help  eliminate  false 
clutter  returns.   Multiple  PRFs  are  used  to  detect  weak 
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targets  that  ordinarily  would  be  masked  by  rain  clutter 
having  the  same  ambiguous  doppler  frequency. 

Because  ground  clutter  is  so  spotty  in  nature,  great 
difficulty  is  experienced  in  calculating  appropriate 
threshold  values  for  detection.   To  solve  this  problem  a 
digital  ground  clutter  map  is  implemented  with  one  word 
for  each  range-azimuth  cell.   All  these  are  stored  on  the 
magnetic  disc  memory.   Sometimes  magnetic-bubble  and  EBAM 
(Electron-Beam-Addressed  Memory)  are  used,  both  of  which 
have  the  advantage  of  being  nonvolatile.   Also,  CCD 
memories  are  used,  in  place  of  disc  systems,  which  are 
suitable  for  storing  the  clutter  map. 

MTD  is  said  to  have  subweather  visibility.   In  fact 
it  can  see  a  small  commercial  type  aircraft  whose  cross- 
section  is  many  dB  below  the  radar  cross-section  of  the 
weather  return.   This  feature  was  not  previously  available. 
Now  it  is  used  to  provide  excellent  tracking  of  aircraft  in 
rain. 

A  second  feature  of  the  filter  bank  approach  is  the 
facility  with  which  proper  thresholds  can  be  established 
taking  into  consideration  the  presence  of  rain.   The 
detection  threshold  is  established  by  summing  the  detected 
output  in  range  cells  of  interest. 

The  MTD  has  the  advantage  of  being  able  to  detect  and 
track  zero  doppler  aircraft  targets  while  ignoring  bird 
and  insect  targets.   At  zero  doppler,  aircraft  targets  are 
broad-side  and  hence  have  much  larger  cross  sections  than 
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do  birds  or  insects.   Therefore,  a  higher  threshold  can 
be  used  in  a  zero  doppler  filter  to  eliminate  bird  echoes 
but  still  detect  aircraft  targets. 

Tracking  is  automatically  initiated  and  the  tracks  are 
continuous  despite  ground  and  weather  clutter  and  birds. 
The  radar  now  can  be  sited  freely  without  consideration  to 
ground  clutter  limitations. 

Both  MTI  and  pulse  doppler  processing  have  the  disadvan- 
tage in  regard  to  vulnerability  to  jamming.   For  both 
techniques  many  more  than  one  pulse  detection  is  required 
to  see  a  signal  imbedded  in  clutter. 

Finally,  no  tuning  or  adjustment  of  any  kind  is  needed 
for  this  combination.   It  is  a  reasonably  economical 
solution  to  the  radar  automation  problem. 

C.   FFT  APPLICATION 

Within  the  past  several  years,  digital  processing  of 
radar  signals  has  become  quite  popular.   Recent  advances  in 
A/D  converters,  integrated  technology,  and  the  emergence  of 
the  Fast  Fourier  Algorithm  have  made  practical  the  digital 
filtering  of  signals. 

The  application  of  the  FFT  algorithm  starts  from  the 
first  known  technique  of  Convolution  by  the  Discrete 
Fourier  Transform.   It  is  shown  in  Figure  16,  where  a 
block  diagram  explains  the  process. 
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Figure   16.      Convolution  by  discrete  Fourier  Transform 
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It  is  simple  in  concept  and  consists  of  taking  an  input 
sequence  of  samples  from  the  received  signals.   Then  a 
Fourier  Transform  is  performed  to  yield  samples  in  the 
frequency  domain,  multiplying  these  by  a  stored  filter 
function  and  performing  an  inverse  transform,  yields  the 
time  domain  response. 

Thus,  any  filter  function  can  be  implemented  by  storing 
it  in  memory.   So,  knowing  the  clutter  spectrum  one  can, 
in  principle,  filter  the  clutter  by  implementing  the 
appropriate  function. 

The  FFT  algorithm  is  shown  in  the  pipeline  of  Figure  17, 
This  technique  is  used  to  digitally  realize  a  bank  of 
contiguous  filters  and  thus  produce  a  spectrum  analyzer  or 
doppler  filter  bank. 

Advantages  of  the  system  include  stability,  reproduci- 
bility, and  flexibility.   The  amount  of  time  for  this 
processing  is  a  parameter  of  considerable  interest  because 
real  time  solutions  are  commonly  required.   To  reduce  time, 
high-speed  A/D  converters  are  being  developed  and  employed 
in  modern  systems  . 

Following  the  discussion  of  the  several  techniques 
described  in  Chapters  II,  III  and  IV,  a  tabulation  in 
matrix  form  comparing  filter  characteristics  with  clutter 
characteristics  has  been  prepared  and  is  shown  in  Table  I. 
The  classifications  of  techniques  are  qualitatively  identi- 
fied as  Superior,  Intermediate  and  Inferior,  in  terms  of 
the  amount  of  clutter  suppression  or  improvement  factor  in 
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dBs  and  also  the  range  at  which  each  process  is  reasonably 
sensitive  where  data  information  is  available. 

The  purpose  of  this  comparison  is  to  help  in  making  a 
decision  about  the  system  to  be  selected  for  a  give  use. 
It  is  noted  that  the  decision  will  be  based  on  the  expected 
requirements  from  the  system  the  area  where  it  is  to  be 
installed,  if  a  fixed  or  mobile  site,  and  on  other  condi- 
tions which  need  consideration. 

Each  radar  system  is  generally  designed  to  optimize  the 
realization  of  a  specific  set  of  measurement   or  observa- 
tions.  Once  the  design  criterion  has  been  achieved,  usually 
it  becomes  difficult  to  utilize  the  radar  for  other  purposes 
without  extensive  and  costly  modifications. 

The  observations  and  measurements  shown  in  Table  I  were 
obtained  from  a  number  of  independent  observers.   The  data 
was  taken  chiefly  at  X-band  on  appropriately  configured 
Navy  and  air-surveillance  radars. 


51 


<ixi<w<ciij        <w<w     <:ijqrt:tti     <pa     <:w 


^  in       ^  in 
in  ■<3"       in  ^ 


PQ  CQ  PQ  OQ  OQ  CQ 


a 


<X3  -P 

5  3 


b 


9  i 

o  o 

<   00  <  00 

I  I 

o  o 

in  in 


m 


a  OQ  cq  Q 


eel 

CQ  CQ  < 

O  O  O 

CM  CN  in 


OQ  OQ  CQ  CQ  OQ  CQ 


U  U  OQ  U 


9.    S 


2  o  5       o 

<  CQ  CN  OQ  CN         CQ 

O  I    O           I 

in  in  in       in 


CQ 


9 


o 

OQ  CN 

I 

in 


U 


OQ 


to 


CQ  CN         CQ  CN  in      CQ 
in  in   '^ 


€ 


CQ  O  U        <:  in  CQ  U 


VD 


ffjrj.  cq5cqq7^cqquq      cqqoqq      cqq      uq 

o  o 

m  m  ro 

<-{ 

§     §      g 

OQCN  CQ^CCQin  CJ  U  CQQ         U  CQQ  U 

I  I 

in  in  in 

•H  rH  Q 


ro 


w 


I 

u  c^ 


>1 

-p 


4:* 


rH 
U 


u  <^      u  <^ 


I 


o 


I 


s 

w 


u 


.5 


e 


CQ  CQ  OQQ  UQU  OQQCQOCQQ  U 

O  O  <H 

CN  CN  •• 

Q 

o           ^  ^         o 

CQQCQCNCQ  CQ  CQCNUQU           CQQ           U 

I             in  in           I 

O              iH  iH           in 


O 


I 

rH 
U 


u 


I 

-H 

a  c; 

•H   o 

U    fO 
CI*  -P 


u 


CQ  H 


u 


CO 

rH 

0) 
5» 


<:  u 


H 


II     II 

5-1    5-1 

.2.2 

CQ 
< 

vo 

1 

Super 
Infer 

Eh 

H-t 

B 

4-4 

Q) 

^ 

CJ 

> 

52 


V.   AN  APPLICATION  OF  A  DIGITAL  FILTER 

In  the  radars  mentioned  earlier,  use  of  analog  MTI 
systems  give  a  clutter  suppression  of  13.8  to  24.6  dB  for 
rain  clutter  according  to  the  data  of  Table  I.   This  suppres- 
sion is  not  sufficient  and  the  problem  becomes  greater  when 
the  radar  system  is  operating  under  Mediterranean  environ- 
mental conditions,  where  the  fall  and  winter  period  is  full 
of  heavy  rainfalls. 

The  experience  of  the  Greek  Navy  for  certain  clutter 
due  to  sea  and  birds  does  not  concur  with  the  suppression 
values  given  in  Table  I.   A  level  of  18  to  25  dB  is  achieved 
under  perfect  radiation  conditions. 

Since  the  radars  under  consideration  are  in  use  in  the 
Greek  Navy,  a  requirements  exists  for  the  improvement  of 
their  signal  processors  in  order  to  get  an  improved  clutter 
suppression. 

This  improvement  can  be  achieved  by  the  use  of  a  higher 
than  second-order  doppler  filter  according  to  Ref.  12. 
This  then  is  the  motivation  for  the  design  of  a  fourth  order 
recursive  digital  filter. 

For  this  purpose  a  Butterworth  characteristic  band- 
pass filter  is  to  be  used  with  an  effective  bandwidth  of 
12616.63  rad/sec.   Certain  restrictions  are  put  on  the 
filter  design  as  they  are  defined  in  the  design  part  of 
this  thesis. 
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The  analysis  and  design  are  done  according  to  Ref.  10, 
The  concept  of  predestorted  or  prewarped  analog  transfer 
function  is  used.   To  get  a  one-to-one  mapping  in  the 
z-domain  approximately  the  same  as  that  of  the  s-domain, 
the  well-known  bilinear  z-transformation  is  used. 

To  detect  moving  targets,  the  system  uses  the  doppler 
shift  in  frequency,  which  is  given  by  the  relation 


2V 

^d    =    -ir  (1-1) 


The  wavelength,  X,  of  the  transmitted  pulse  might  also  be 
used  as  a  parameter.   Since  the  design  of  the  above  filter 
refers  to  specific  systems,  where  A  is  fixed  by  design 
parameter,  this  option  is  not  open.   These  fixed  values  of 
X  for  the  two  radars  are  calculated  by  using  the  relation 


A   =  f-  (5.1) 

o 


where 


g 

c   =   3x10   m/sec,  velocity  of  light 

f   =   operating  frequency  of  the  system  given  in  Hz 


For  f   =  1300  MHz  (AN/UPS-1)  A  =  .23  m; 

and  for  f   =  500  MHz  (AN/SPS-40BC)  A  =  .6m. 
o 

V   stands  for  relative  velocity  of  the  target, 
o 

Usual  types  of  moving  targets,  in  the  area  of  interest, 
are  given  in  Table  II   with  their  cruise  and  maximum 
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Type 

Cruise 

Helicopters  Generally 

80-130 

C-47 

150 

C-130 

200 

ALBATROS 

200 

Commercial  Air 

Lines 

300-400 

T-2E 

250 

A-7H 

300 

F-84F 

350 

F-5 

400-550 

F-102 

500 

F-104 

400-550 

F-ICG 

500 

F-4 

500 

Maximiim 

180 

250 

300 

350 

500 

600 

600 

600 

700 

800 
1000 
1000 
1300 


TABLE  II 


VELOCITIES  OF  SELECTED  TARGETS 
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velocities  in  nautical  miles  per  hour.   This  information 

was  used  in  the  beginning  of  the  filter  design  in  order  to 

evaluate  the  necessary  doppler  frequencies  assuming 

(i)    Radar  stationary 

(ii)   Radar  in  motion  with  a  velocity  varying  from 

5  to  30  nautical  miles  per  hour,  to  correct  the 

effective  bandwidth  of  the  filter. 

Table  III  contains  a  summary  of  the  literature  [Ref.  2, 

Chap.  6;  Ref.  3,  pp.  146,  551-554;  Ref.  19,  p.  17-9]  for 

the  relative  velocities  and  standard  deviations,  a  ,  of 

'   V 

clutter  spectrum  in  meters  per  second.   These  kinds  of 
clutter  are  common  to  the  Mediterranean  environment. 

The  conversion  of  the  standard  deviation  of  the  clutter 
power  spectrum,  a  ,  to  frequency  is  obtained  by  applying 
the  relation 


2a. 
c 


o  =       -r^  (5.2) 


where  A  is  the  wavelength  of  the  transmitted  pulse  in  meters. 
This  relation  is  an  alternative  to  equation  (1.1). 

Computer  programs  2,  3  and  4  were  used  to  plot  the 
amplitude  square  characteristic  versus  frequency,  w,  of 
(1)  analog,  (2)  prewarped,  and  (3)  digital  transfer  functions. 
Figures  19,  20  and  21  show  the  characteristics  of  the  designed 
filter. 
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A.   DIGITAL  FILTER  DESIGN 

The  second  part  of  this  thesis  is  the  design  of  a 
digital  filter  which  is  going  to  reject  clutter  echoes  in 
ranges  indicated  by  Table  III  [Refs.  8,9,10,11,12].   Certain 
numbers  of  doppler  filters  combined  with  range  gates 
constitute  the  video  part  of  an  MTI  radar  which  is  shown 
in  Figure  18  in  block  diagram  configuration. 

A  good  performance  criterion  of  MTI  doppler  filters  is 
the  MTI  improvement  factor.   This  factor  defines  the  improve- 
ment of  Signal  to  Clutter  ratio  which  can  be  achieved  by  a 
filter.   In  the  table  of  Ref.  12,  the  improvement  factor  I 
is  given  as  a  function  of  the  order  of  the  filter. 
Generally  the  higher  the  order  of  the  filter  the  better  the 
improvement  factor.   However,  higher  order  transfer  functions 
imply  complexity,  low  reliability  and  higher  cost. 

The  goal  of  the  present  doppler  filter  design  is  to 
derive  a  transfer  function  which  will  provide  high  clutter 
suppression  and  at  the  same  time  reduce  the  complexity 
caused  by  the  higher  order  filter,  to  a  reasonable  level. 
For  the  above  reasons,  a  fourth  order  digital  filter  will 
be  designed.   The  improvement  factor  will  be  increased  and 
the  filter  complexity  will  be  kept  to  a  minimum 

An  Infinite  Impulse  Response  (IIR)  band-pass  filter  of 
the  fourth  order  will  be  implemented.   As  is  known,  the  IIR 
filter  does  not  need  as  many  delays  as  the  Finite  Impulse 
Response  (FIR)  filter.   Consequently,  the  number  of  weights 
is  considerably  lower. 
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The  transfer  function  for  a  fourth-order  band-pass 
digital  filter  is  given  by  the  following  equation 

R(^)       -      A(^)   -   ^0-^^l^"^-^^2""^^^3^"^-^^4^"^    ,,  ,, 

^^^^  l  +  b,z  -^^-b^z  ^+b,z  '^  +  b,z  ^ 

12       3      4 

This  transfer  function  is  completely  defined  by  the 
values  of  a^,  a^ ,  z^,    a    ,    a.,    h^,    h    ,    b^,    h..      There  are 
two  ways  to  define  the  above  coefficients.   The  first  is 
to  go  directly  to  the  z-domain.   According  to  the  frequency 
response  characteristics  which  will  be  defined  later,  one 
has  to  find  the  exact  position  of  zeros  and  poles  of  H(z). 
Since  zeros  and  poles  are  calculated  the  values  of  a's  and 
b's  occur  after  a  simple  multiplication  of  the  factorized 
polynomials  A(z)  and  B(z). 

The  second  method,  to  be  used  here,  is  to  select  the 
desired  H(s)  transfer  function.   From  this  is  defined  the 
transfer  function  H(z).   It  is  more  difficult  than  the  first 
method  but  has  the  advantage  of  permitting  design  of  the 
filter  in  both  the  s  and  z  domains. 

The  concept  of  predistortion  applied  to  the  transfer 
function  H(s)  is  introduced.   To  transform  H(s)  to  H(z) 
the  z- trans formation 


-sT 
z   =   e    ^  (5.4) 


is  normally  used,  where 


T   =  sampling  time, 
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Because  of  the  cumbersome  calculations  resulting  in  the 
use  of  (5.4),  the  simpler  bilinear  z-transform 


s 


is  used.   This  is  an  approximation  to  (5.4) 

The  use  of  bilinear  z-transform  distorts  the  frequency 
response  characteristics  of  H(s) I   .  .   For  this  reason 
H(s)  is  predistorted  in  such  a  way  that,  after  the  bilinear 
transformation,  H(z)  will  have  the  desired  frequency  response 
characteristics . 

Starting  in  the  s-domain  a  Butterworth  fourth  order 
band-pass  filter  is  selected  because  of  its  flat  response 
characteristics.   The  transfer  function  of  such  a  filter  is 
given  by  [Ref.  11]: 

2  2   2 

H(s)   = 3 (5.6) 

s   +a  —  s   +(2  +  -2)a)^S+a  -Q-  s  +  oJq 

V  Q  « 

where    k  =  gain 

u)   =  27Tf   =  center  circular  frequency 
o      o 

Q  =  quality  factor  which  is  given  by  Q  =  u  /BW  and 

BW  =  3  dB  bandwidth 

a  =  1.41421 

(constants  given  in  Tables  of  Refs. 
b  =  1.0       7  and  11) 
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If  BW  and  oj^  are  given,  the  frequencies  at  the  3  dB 
points,  f^   and  f^,    can  be  defined: 


^^2  "  "l  ^      ^^^2  "  ^""^1  "   ^^  ^^-"^^ 


2 

Wj^CO^    =    27Tf^-2Trf2    =    0)^  (5.8) 


The  following  restrictions  are  imposed  on  the  filter 
design  using  charts. 

(i)   The  maximum  doppler  frequency  for  clutter  (angels, 
birds,  insects)  is  100  Hz.   This  corresponds  to  a  velocity 
of  58.31  nm/h  for  X  =  .6  m  and  22.35  nm/h  for  X    =    .23  m. 
The  desired  clutter  suppression  at  100  Hz  is  30  dB. 

(ii)   The  filter  is  to  be  used  in  air  search  radar  and 
must  detect  targets  with  velocities  of  881.7  nm/h  for 
X  =  .6  m  and  340  nm/h  for  A  =  .23  m.   This  sets  the  upper 
limit  of  frequency  at  9500  rad/sec  or  1572  Hz. 

(iii)   Slow  moving  targets  such  as  helicopters  must  be 
detected  with  a  suppression  of  no  more  than  3  dB.   This  sets 
the  cut-off  frequency  at  this  end  of  the  spectrum  at 
400  -  600  Hz. 

With  these  restrictions  the  following  values  were 
selected. 


0),  =  2829.44  4ad/sec,   f   =  445.63  Hz 


(jO 


=  15500  rad/sec,  f^   =  2463.7  Hz. 
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Then,  applying  basic  relationships 

0)^  =  6622.47  rad/sec 

Q   =   .525 

BW  =  12616.63  rad/sec. 

With  the  constraints  above  and  these  values  the  function 
is  prewarped  and  the  bilinear  z-transform  applied.   The 
final  filter  parameters  were  obtained  after  successive 

applications . 

1      1 2 
In  Figure  20  the  prewarped  form  of  |H(jcj)  |   vs .  u  is 

plotted.   The  prewarping  was  made  on  the  center  frequency 

uj   with  the  exact  Nyquist  sampling  time,  T   =  .2479445x10 

sec. 

In  Figure  21  the  final  filter  power  spectrum  is  given 

after  the  application  of  bilinear  z-transform.   Comparing 

this  curve  with  the  analog  plot  in  Figure  19,  it  may  be 

seen  that  a  considerable  distortion  exists  for  frequencies 

(jj  >  CO  .   This  is  reasonable  as  a  result  of  the  application 
o 

of  the  transformation 


2  "d'^ 

0)    =   =—  tan  — X —  (5.9) 

a     T        2 

s 


where    w   =  analog  circular  frequency,  rad/sec 
a 

CD   =  digital  circular  frequency,  rad/sec. 
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For  small  values  of  w^,  to^  =  im^,    but  for  high  values  of 

'^a'  '^a  ^  '^D*   "^^^s  is  implied  by  the  nonlinearity  of 
equation  (5.9) . 

The  final  filter  design  parameters  are  obtained  from 
Figure  21  and  are  as  follows: 

0)2  =  9900  rad/sec 

f^  =  1575.63  Hz 

w,  =  3465  rad/sec 

f^  =  551.47  Hz 

Aio  =  6435  rad/sec 

Af  =  1024.16  Hz. 

More  details  of  the  analysis  are  given  in  Appendix  A. 
Clutter  suppression  of  34.2  dB  was  calculated  at  the 
frequency  of  100  Hz.   This  is  better  than  that  of  the  analog 
delay  line  which,  in  the  existing  MTI  radar,  is  between 
25  and  30  dB. 
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VI.   SUMMARY  AND  CONCLUSIONS 

This  thesis  was  motivated  by  a  requirement  to  find  an 
improved  anticlutter  technique  for  MTI  radars  presently  used 
in  the  Greek  Navy.   Initially,  a  complete  survey  of  existing 
methods  of  clutter  suppression  was  undertaken.   Table  I 
summarizes  the  results  of  this  survey  showing  qualitatively 
clutter  improvement  against  various  methods  and  types  of 
clutter. 

The  Fast  Fourier  Transform  (FFT)  appears  to  be  the 
optimum  selection.   However,  it  does  not  lend  itself  easily 
to  introduction  into  the  older  radars  still  in  use.   There 
is  no  doubt  FFT  is  used  in  modern  radar  systems  and  will  be 
extended  in  the  future,  where  cost  and  trained  personnel  are 
no  obstacles. 

For  the  application  under  consideration  a  method  of 
digital  signal  processing  was  selected.   It  has  the  advan- 
tages of  being  practical,  reliable  and  flexible.   In 
addition,  due  to  advances  in  LSI  technology  the  digital 
filter  may  be  contained  in  a  small  package  of  light  weight 
with  good  performance  characteristics  and  can  be  constructed 
at  reasonable  cost. 

The  design  idea  used  here  was  introduced  to  the  author 
for  the  first  time  in  a  digital  filter  course  at  the  Naval 
Postgraduate  School. 
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The  concept  of  the  prewarped  transfer  function  was 
employed.   This  approach  permits  one  to  carry  out  his  design 
in  the  familiar  s-domain.   The  design  constraints  were  taken 
from  Table  III.   This  table  was  assembled  from  a  number  of 
sources  giving  clutter  characteristics  for  an  unusual 
variety  of  clutter  types.   Since  these  —  birds,  insects, 
angels  and  weather  (squalls)  —  are  typical  of  the  ones 
found  in  the  Mediterranean  environment  against  which  the 
present  radars  do  not  give  satisfactory  suppression. 

The  transfer  function  H(s)  was  chosen  from  a  table  of 
filters  of  the  Butterworth  type.   It  is  a  fourth  order 
filter.   Knowing  H(s),  predistortion  was  applied  in  order 
to  compensate  for  the  distortion  introduced  by  the  applica- 
tion of  the  blinear  z-transform.   This  is  a  powerful  digital 
filter  design  tool  since  it  permits  the  designer  to  control 
at  each  step  the  values  of  the  filter  parameters. 

This  thesis  actually  carried  out  the  preliminary  design 
of  the  selected  filter.   The  results  were  then  plotted 
using  Continuous  System  Modeling  Program.   The  filter 
parameters  used  in  the  program  were  those  derived  from  the 
clutter  characteristics  surveyed  as  part  of  the  thesis 
study.   At  the  100  Hz  cut-off  frequency  this  filter  gives 
4.2  to  9.2  dB  improvement  in  clutter  rejection.   This  is 
a  significant  improvement  over  the  existing  DLC  MTI  systems 
and  is  considered  to  have  met  the  clutter  requirements  for 
ships  operating  in  the  Mediterranean  environment. 
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The  adaptive  nature  of  the  design  proposed  for  the 
clutter  problem  suggests  that  it  may  be  further  utilized 
for  other  applications  by  interested  students  in  the  future, 
Additionally,  digital  filters  may  be  actually  constructed 
using  the  design  program  developed  here  and  tested  for 
correlation  with  the  design  parameters. 


71 


J" 


0) 

c;^ 

c 

(0 

u 

• 

> 

S-l 

0 

-p 

u 

t-i 

d 

y-l 

5 

4J 

2 

c 

o 

OJ 

o 

e 

A 

OJ 

Pm 

> 

0 

« 

u 

u 

a 

Eh 

e 

D 

•H 

Pu 

s 

Sh 

o 

o 

u 

4J 

JJ 

0 

iH 

0 

IM 

0 

■p 

0 

Oi 


2. 


.-4 

J^ 

• 

—t 

-H 

II 

>-■ 

» 

_J 

—. 

u-^ 

•-t- 

O 

"■"' 

O 

• 

a. 

M 

r 

» 

^ 

X      - 

< 

• 

• 

'J5 

a: 

■^ 

iT* 

>— « 

o 

• 

-^ 

il 

;r\ 

II 

i. 

Cl. 

J^ 

1— 

•» 

CL 

• 

-; 
-< 

O 

-J> 

•-• 

'-3 

II 

^ 

• 

-w 

» 

r- 

»— » 

• 

■■>J 

• 

J— 

-J 

tn 

■ff 

r-4 

3 

UJ 

'JO 

» 

tl 

t-^ 

Q 

• 

QC 

-J 

"S 

o 

» 

~« 

UJ 

«k 

s: 

JO 

* 

'-3 

-^V- 

• 

>r 

J— 

oo 

sr 

» 

rj 

O 

•  J 

^-•i 

-^ 

r^ 

Oq^ 

H- 

»«» 

• 

» 

o>- 

^ 

» 

--rf 

• 

ox 

> 

» 

•* 

O 

!NJ 

/O 

>-. 

^j  — 

M 

•  'J 

• 

•  «M  Nl  'O  — • 

• ',» 

y) 

^•^'V- 

''_ 

>    !l 

O-^ 

•D 

II 

0   1 

■s. 

u  %. 

X 

-•■A. 

■  D 

< 

<>  <» 

—4 

T.  -  - 

i. 

^ 

r— J     -^» 

^^ 

».— •  ••- 

/* 

■"O 

z: 

I 

II            <l 

H 

t-   -" 

^ 

"^     • 

^ 

-<       %j 

^ 

ii  >-« 

-/; 

•  ct 

I— 

f— 

o--i 

-  J 

OlU. 

OL 

2 

•X; 

O 

,>~_ 

< 

:x: 

X  ^.:  'X         CL 

«/> 

-^ 

^- 

; — 

»-!—       '■^'^ 

a^ 

vi 

_< 

JXO       T\- 

sA 

CL 

V 

■< 

•LO      --ir 

72 


CO 


tP 

o 

rH 

(0 

c 

cd 

U-l 

0 

fvi 

o 

5" 

•H 

2 

JJ 

2 

U] 

o 

•H 

o 

u 

a 

<D 

CU 

4J 

u 

« 

fO 

w 

u 

E-t 

(0 

D 

x: 

04 

o 

2 

O 

OJ 

U 

^3 

P 

+1 

•H 

iH 

CU 

(0 

+> 

o 

04 


It* 

iA 

<*» 

V- 

o 

a 

t- 

3 

C 

a: 

o 

H 

<r 

-J 

v> 

;i£ 

< 

o: 

1- 

• 
0 

• 

f*i 

0 

0 

2. 

0 

(T 

•-< 

c 

H 

> 

0 

-J 

a. 

^_ 

in 

0. 

1— « 

<m. 

» 

» 

1— • 

• 

c 

• 

i 

rs»o 

fO 

X 

« 

wc 

H 

>— 

UL 

rHO 

1- 

CJ 

Ui 

■»•  J^ 

^ 

►-< 

(^ 

U-O^ 

U4 

^ 

u. 

♦  a\ 

LJ 

# 

u. 

U-l"^ 

u\u\ 

» 

# 

*  'vO 

t    • 

• 

p- 

O 

OiSt 

»  «■ 

0 

11 

^ 

*.Cj 

-ft  » 

0^ 

H 

IT  ♦ 

.•^^-> 

tl 

0 

U. 

<— •  '. 

<TJ<-J 

LU 

»-^ 

•  J 

* 

u- 

Q 

LU 

ft 

'JU* 

32  LJ 

^ 

•> 

• 

^tu. 

COO 

CJ 

1- 

i 

«^ 

>J-ii 

<^'<^ 

»-t   • 

o 

LU 

tn 

rj:0 

+  -♦" 

:x  • 

_J 

s> 

i»- 

1  i^ 

^^.^t 

->o 

Ol 

• 

(JV 

u.    • 

<IOfNX<J 

>■ 

> 

r- 

*  r>j 

<o>-*o 

X 

</) 

.-» 

u.>r 

*  «  ^.r^tn 

a* 

*  CO 

<o 

r-4:s>r 

Ui 

</} 

Lf\ 

u-r«- 

<C>"'n 

0 

c 

1 

1*  t— ( 

II 

U     (1     II 

>-<  n  " 

3^>-<-n 

•  J. 

iTfcoc: 

{-■^CM 

n-j^:ji3  » 

< 

■■coc. 

l1>->XXf- 

at^U. 

1- 
2: 
a 

u. 

u. 

U 

J. 
1- 

a: 

CK-O 

ij 

j- 

4CCL 

0 

•A 

3: 

* — 

!-»»- 

V-l 

iJUnCO 

;— 

XQ-C 

a 
oc 
zs- 


73 


o 


(0 


QC 
O 


TJ 

(U 

a 

}-i 

to 

5 

Q) 

S-i 

a 

H-l 

^*n 

0 

en 

o 

5" 

■H 

3 

4J 

2 

cn 

u 

-H 

o 

>-i 

« 

0) 

CU 

4J 

u 

DS 

(0 

U 

M 

Eh 

(0 

D 

£ 

CU 

o 

S 

o 

0) 

u 

'O 

3 

+J 

•H 

rH 

CU 

6 

Id 

+) 

0 

04 


< 
a. 

I- 


<u 


o 

a- 

M 


Uu 

i? 

• 

(-J 

in 

o. 

•-4 

U^ 

a 

»-* 

ou. 

» 

*• 

•-H 

c# 

• 

• 

o** 

c 

r- 

mlr^ 

»-4 

n 

2. 

>oo 

m 

i 

< 

" 

m-*- 

M 

b- 

a. 

^gj 

»— 

u 

i^ 

■♦■  • 

-J 

►— • 

Ul 

U.!"* 

^ 

J. 

K 

H-VJU 

u 

» 

Lk. 

u.>y^ 

UIUA 

v> 

• 

«  p- 

•     • 

« 

r^ 

o 

U30 

tr  * 

o 

II 

z 

•  >n 

*  « 

0^ 

t~ 

1-4 

'jp» 

i^>M» 

It 

X. 

_J 

U- 

^.-« 

— «« 

^ 

o 

lU 

* 

yJJW 

Cf2«-) 

UJ 

l-< 

o 

u. 

vC    -♦• 

CCC_) 

Q 

UJ 

(_; 

■IT 

WIX. 

*  * 

r- 

^. 

21 

l-< 

(M* 

CDC_) 

-J 

» 

• 

CNJV-U 

OJU 

a 

1- 

5: 

■O 

C-iii 

-^•^ 

—I  • 

c 

f  ■ 

•c 

-Til- 

+  •♦■ 

:i  • 

-J 

>- 

L-C 

1  «■ 

—.-• 

->o 

c 

V, 

0* 

u.    • 

<i.O«V'J-U^ 

>- 

> 

rs) 

»  l-O 

<o 

>-*  03 

X 

t/J 

(M 

turn 

»*  -sr-IJU 

fM 

»m 

<o 

r-ITE,-! 

yj 

V7 

r- 

U-rr, 

<IO>-~)  11 

o 

_; 

rsi 

« a  (\j 

-WW. 

II    i».>t 

^< 

o 

1 

C'J-  t 

2.—^ 

i-t  !I  — .U-LL-JO. 

3 

II 

II   II    II 

h-*  II    II  JC^f->/7 

•"X 

2 

<cL^.^a 

•—r-l(\i-i-}^^S.     • 

>— « 

<i 

cC^C 

l(  >->-X-U>-.ci:-5u_ 

}- 

U- 

U.       X 

2 

Q 

>— 

a 

a:c^-I^ 

tj 

_  J-  ^  cr         a 

M* 

SJ-> 

->-     oe 

74 


^ 
A 


J) 


ta 


04 


< 

_J 

a 

I/- 

0 

i* 

<\i 

< 

u^ 

Q. 

*«.»'• 

•"» 

»- 

it:— 

• 

r> 

*  — • 

*  "^ 

u 

2: 

•  u_ 

•  u. 

tr 

r— 1 

C'* 

n* 

II 

> 

«t 
li' 
0 

•-^ 

?f  t- 

^;  <^ 

a 

.- 

•4-r- 

sr« 

•> 

•> 

>-- 

•x>r 

ODsj- 

• 

• 

•    ceo 

Coo 

c 

r- 

»-<    ^-»r 

r-^ 

i-t 

11 

z. 

•f     fvjm 

ovjrn 

m 

X 

< 

->—  .o 

•0 

H 

>- 

ac 

u-u-in  • 

ay   • 

>- 

0 

0 

«•  fl^  — "-^ 

*^r-« 

^ 

•-« 

LJ 

iit^  +  '~' 

"♦"- 

— U. 

JC 

U. 

*  f»  — f 

'"^ 

-•u 

♦ 

U. 

•  »— '^^ 

•»»i« 

u^u^ 

-*  • 

• 

(NJCMU-'^ 

LL  — 

1   • 

->  • 

r- 

V^ 

**»•»■»  a. 

■*ll. 

♦  ir 

xo 

M 

iZ 

(>0 -?:*:* 

:<:♦ 

*  ■* 

•W»7> 

t- 

•~t 

O'-*  it 

♦  it 

,«>••-> 

r-«  II 

X 

_j 

o</i  •* 

•  « 

*i»— 

:x-i 

0 

U.' 

•R-  *  >r  • 

>r  • 

cuo 

-?UJ 

»-< 

Q 

«     twoj 

•~-fvj 

OJO 

XC 

x 

LJ 

04v>,j(/.*' 

^'-' 

•*  "^ 

— 'r- 

X 

2. 

a^  I     1  U*^ 

«— ^ 

ojlj 

*I3 

» 

>5--.-^OU 

V5W 

OJCJ 

trvo 

r- 

2: 

^r'j.!x«-t_; 

*  (/I 

^•'-mf 

r*  • 

0 

C  «   -i  (J-  » 

(^» 

+  -♦• 

tvj  • 

-J 

y- 

r-i<:^<;Nl- 

O^t 

o**^ 

u>o 

a. 

wO 

sT  *  5t  =0  -^ 

aOr-4 

^OCNJI—UA 

P" 

> 

fN)    t    •lO-O 

if\>C 

<%<^ 

>>rco 

X 

v^ 

OsT^f-^"^ 

_jU-cO^ 

it* 

-N.LT'CO 

•;^--.^   .,-n 

_j  •ro^ 

<o 

r-»     ••-< 

Ul 

o") 

.-(/)p"-f^    » 

+  ir.  • 

4- 

<o>tn  II 

0 

3 

0 

"•Q:  <irirv.'a.r<^r\i 

r> 

v.'J^^'"' 

«l  U^SX 

3< 

•1  .^ijT.^.^ 

CI.  *^*-* 

r-> 

21—''-' 

t-i    11   1— 

u-u-na. 

3 

jC  n  11  II  w 

II  It  II 

It 

1-^  tl    II 

jCiJ- 

vj 

•X 

z 

<j-j2a.-)oO  jc 

1^ 

.^.^(Nj-)-);^ 

^3   * 

»— 1 

;— <rCQ.  _luO  ^ 

r-\ 

11  >->XX^- 

li 

->x 

z 

< 

u. 

LL 

C 

X 

U 

1- 
0 

■•-ii.    <•-       21— 

'  -  T  a  c     '^  «/^ 

75 


APPENDIX  A 
Analysis  of  Digital  Filter  Design 

This  appendix  covers  the  analysis  in  detail  of  the 
digital  filter.   The  frequency  domain  transfer  function  is 
given  in  Ref.  11,  pge  10-4. 

H(s)   =  a)^^bksVQ^ 

o'^  ^   ^o  o3  ^  ,o  ^  b  ,   2  2^    4        ^^'^^ 

S   +   -r-  S   +   (2+—;^)  CO   S   +  OJ 

Q  -2   o       o 

where 

k  =  1,  gain 

a  =  1.41421 

(constants  from  Table  1,  p.  10-4  of  Ref.  11) 
b  =  1.0 

u   =  6622.47 

(defined 
Q  =  .525     in  the 

analysis) 
BW  =  12616.63 

Applying  these  values  to  (A.l) 

H(s)  =  1.59x10^  s^ 

s"^  +  1.784xl0'*s^  +  2.468xl0^s^  +  7 .  824xl0-^-'-s  +  1.923xl0-'-^ 

(A. 2) 
Applying  s  ->  ju     (A.2) 

^,-^    =  -1.591x10^  u?  +   jO 

(0)^^-2.468x10^00^+1.923x10^^)  +  j  (-1.784xl0^co^+7.823xlO-'--'-oo) 
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This  function  was  programmed  in  the  computer  and  the 
results  plotted  in  Figure  19.   Associated  program  was  Program 
Number  2 . 

According  to  the  prewarping  concept,  the  center  circular 
frequency  oj^  is  prewarped.   Since  BW  =  12616.63  rad/sec, 
then  Af  =  2008  Hz,  and 


Tg   =  ~j     =   .2479445x10"-^  sec  (A. 4) 


The  distortion  factor  oj '  is  given  by  the  relation 


u)'   =•  -£-  tan  ^  (A. 5) 

s       ^ 


But  predistortion  means 


Scu 

—?-  (A. 6) 


Applying  (A. 6)  to  (A.l)  and  substituting  the  known  values 


2  2 

li'  (s)       =      (co'/Q)   s 

s"^  +  1.41421(^')s^  +  (2  +  i  (oo')  V  +  ^■^^f^'(a3')-^s  +  (o)')^ 
^  Q  ^ 

(A. 7) 


Substituting  (A. 4)  to  (A. 5)  becomes 


=  8662.18  (A. 8) 
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Since  Q  =  .525,  (A. 7)  becomes 


H'(s)   =  —, 2.12^Q     s 


s^  +  2.33x10^3^  +  4. 2x10  V  +  1.7xl0-^^s  +  5.63x10-^^ 

(A. 9) 

Applying  s  ->  jco   (A.  9)  implies 

„,/.  N   _ -2.72x10^0)^  +  10 

^^'^'       ~     ~A 8~2 15 4"1 12~ 

(oj  -4.2x10  to  +5.63x10^)  +  j  (-2.33x10  co  +1.7x10  co) 

(A. 10) 

This  function  was  programmed  to  the  computer  and  results  in 
Figure  20.   Associated  program  was  Program  Number  3. 

Applying  the  bilinear  z-transform,  which  is  given  by 
the  relation 


s 


to    (A. 9)    becomes 


H(z) 


z     -  2z     +  1 


35.95869z'^  +  5.27084z^  +  2.3814z^  +  1.03484z  +  6.19269 


(A. 12) 

Finally,  this  function  was  programmed  in  the  compuer  and 
the  results  are  plotted  in  Figure  21.   Associated  program 
was  Program  Number  4 . 
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